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Abstract 

The use of energy flow is investigated as a diagnostic tool for determining the color 
SU(3) representation of new resonances. It is found that the pattern of soft gluon 
radiation into a rapidity gap depends on color flow in the hard scattering, and reflects 
the gauge content of new physics. The massive soft anomalous dimension matrix for 
rapidity gap events is introduced for describing soft gluon emission analytically in 
heavy quark pair production. A gap fraction is used for quantifying the amount of soft 
radiation into the gap region. In general, the results illustrate that radiation is greater 
for a singlet resonance than for an octet. Especially, it is found that the quantitative 
difference is quite distinguishable for spin-1 resonances, depending on the gauge content 
in the new sector. 



1 Introduction 

QCD has been successfully tested in many hadronic processes at high energy. Current 
efforts in perturbative QCD have moved toward precision, to help gain accurate predictions 
for Standard Model processes in collider experiments as backgrounds to new physics. Precise 
prediction is a key to find new physics, because QCD events should be distinguished from 
signals in order to claim the discovery of new phenomena at colliders. In this paper, we show 
how to use factorized cross sections as a tool for analyzing the gauge content of new heavy 
particles. 

Factorization is a statement of the quantum mechanical incoherence of short and long 
distance physics, and plays an important role in the use of perturbative QCD. Perturbative 
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calculations of factorized partonic cross sections depend on the separation of regions of mo- 
mentum space. The three basic regions describe off-shell "hard" partons at short distances, 
energetic, on-shell "collinear" partons near the light cone, and long-wavelength, "soft" par- 
tons (see Ref. [1] for reviews and references therein). Once we obtain a factorized form of 
the cross section for a certain process, resummation can be achieved from the independence 
of physical quantities from factorization scales, which leads to evolution equations [2]. The 
evolution equations can be solved to resum large logarithmic corrections in terms of pertur- 
bative anomalous dimensions. Especially, the soft anomalous dimension matrix introduced 
in Refs. [3,4], which is the main interest in this paper, allows us to describe soft gluon 
radiation systematically in hadronic collisions. 

In the mid 1990's, dijet rapidity gap events with anomalously low radiation in a wide 
interjet rapidity region for electron-proton collision were observed at HERA [5-7] and for 
proton-antiproton scattering at the Tevatron [8]. Such events were predicted from the ex- 
change of two or more gluons in a color singlet configuration, which avoids color recombi- 
nation between jets [9]. A quantitative explanation of these events in terms of factorizable 
cross sections was studied in Ref. [10], which discussed the dependence of the dijet cross 
section on energy flow, Q c , into the central region, in terms of a soft function that is a ma- 
trix in the space of the possible color exchange at the hard scattering. A "gap" here refers 
to low energy flow Q c . Since then, there have been various studies of gap events from this 
point of view. The soft anomalous dimension matrix, denoted by Ts below, for rapidity gap 
events, was calculated in Ref. [10], and used to resum the leading logarithmic contributions 
in Refs. [11,12]. The soft anomalous dimension matrix for three jet rapidity gap processes 
was calculated in Refs. [13-15]. 

Cross sections computed from the soft anomalous dimension matrix organize "global" 
logarithms, following the terminology introduced in Ref. [16], and as such do not organize all 
logarithms of the gap energy. Non-global logarithmic contributions, mis-cancellations from 
secondary emission, discovered by Dasgupta and Salam [16], were resummed in the large 
N c approximation [17]. More recently, Forshaw, Kyrieleis and Seymour have discovered new 
logarithmic contributions which they name "super-leading logarithms", which do not seem 
possible to exponentiate by current techniques [18-20]. The contributions of these non-global 
logarithms, however, are generally not dominant [21,22]. 

In this paper, we study the pattern of soft gluon radiation into a rapidity gap for the 
production of heavy particles in hadron-hadron scattering at leading order in the soft matrix, 



2 



Ts, which resums global logarithms. The pattern depends on color flow in the hard scattering, 
and reflects the gauge content of new particles that participate in the process, for example, 
as s-channel resonances. 

The results of this paper illustrate the use of energy flow, treated by perturbative QCD 
and factorization, as a diagnostic tool for determining the color SU(3) representation of new 
resonances, and the use of resummation for more precise cross sections. Previous work on 
soft gluon radiation for the study of new physics includes Ref. [23] for Higgs boson production 
in association with two hard jets, and Refs. [24, 25] for gluino and squark pair production, 
which describe resummed cross sections to provide more accurate predictions for signals. 

There are many proposals for extensions of the Standard Model. These could provide 
solutions to the hierarchy problem, dark matter explanation, the origin of mass, etc. The 
most common feature of the models is that they include new heavy particles transforming 
under definite gauge groups. In many scenarios, we expect that they are created with O(TeV) 
mass via s-channel processes in hadron-hadron scattering. Possible processes with various 
spins and gauge groups of new heavy particles are studied and summarized in Ref. [26], 
which considers new sectors that decay into a top quark pair. There have been studies for 
identifying signals, including hadronic or leptonic decays from gauge Kaluza-Klein (KK) 
resonance states and from QCD background [27-34]. These studies show that the ability 
to distinguish signals especially from hadronic decays, which have much larger branching 
ratio than leptonic decays, depends on both how strongly the new sector is coupled to the 
Standard Model and to what extent signals and backgrounds are separated [35-42]. Also, 
various methods for spin measurements of the intermediate new particles have been suggested 
and discussed (see Ref. [43] for reviews and references therein). 

Our study begins from this point. Once we succeed in distinguishing a resonance from 
backgrounds, a remaining task is to determine the SU(3) color gauge content of the resonance 
particle, denoted as V. For simplicity we consider below two possibilities; color-singlet (Z') 
and color-octet (G). Since we consider QCD radiation with total energy much less than 
resonance energy, our considerations do not depend on the details of the model. In this 
paper, we assume that the decay width is large enough that a hard function is an effective 
vertex. We will, however, discuss the consequences of a narrow resonance at the end of 
Sees. [3] and [6j 

In many models there are gauge interactions whose coupling with a top quark pair is 
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enhanced in particular. These include models with a KK excited graviton, weak and strong 
gauge bosons in extra dimensions [44,45], as well as heavy spin-0 particles in MSSM and 
two-higgs-doublet model (2HDM) [46,47] and spin-1 coloron and axigluon [48]. Such new par- 
ticles V of mass M and spin (sp) could show up as resonances in pp(p) — > V'(M, sp) — > QQ 
at the Tevatron and the LHC. We may expect larger branching ratios for heavy quark pair 
production than for other channels such as massless dijets or dileptons, due to the small cou- 
plings of V to light particles. This motivates the calculation of an analytic form describing 
soft gluon emission in heavy quark pair production. At the level of global logarithms, this 
is determined by the soft anomalous dimension matrices for rapidity gap processes involving 
massive particles. The "massive" soft anomalous dimension matrix for heavy pair production 
was studied in Refs. [49-57] for threshold resummation. In this paper, we extend this work 
to introduce the massive soft anomalous dimension matrix for rapidity gap cross sections, 
which will be used as a tool for our analysis. 

Reference [10] studied rapidity gaps for a process in which t-channel exchange was dom- 
inant. The result indicates that, in the limit of a very large interjet region, the color singlet 
component dominates, with more radiation into the gap when the exchange is a color-octet. 
In this study, however, we might expect a different pattern of soft gluon radiation due to new 
heavy particle resonances in the s-channel rather than t-channel. The quantity of interest 
is the gap fraction, the ratio of the number of events for heavy quark pair production with 
a gap energy up to some value Qo to the total number of events for the process. We will 
find a quantitative difference that is quite distinguishable and can be used for the study 
of the gauge content in the new sector. This difference, especially caused by the SU(3) 
color content, is related to "drag" effects in e + e~ three-jet (qqg) events, compared to qq^y 
events [58-60]. The color dipole configurations explain a surplus of radiation in the qg and 
qg interjet regions for three-jet events, while radiation is enhanced in the qq region for qqj 
events. 

In this study, we do not give a detailed study of non-global [16] or super-leading loga- 
rithms [18], arising when gluons radiated outside the gap emit back into the gap. The latter 
are sometimes called secondary emissions. We can in principle introduce correlations be- 
tween energy flow and event shapes that suppress such secondary emissions naturally. This 
technique was introduced in Ref. [61] for rapidity gap e + e _ dijet events. Such correlations 
control secondary radiation, by suppressing states with radiation outside the gap at interme- 
diate energy scales, leading to perturbatively computable measures insensitive to non-global 
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logarithms. In our study, the rapidity gap we choose does not allow radiation of high trans- 
verse momentum outside the gap by definition, which suppresses the non-global logarithms 
from final-state radiation. On the other hand, another issue arises for hadronic collisions, 
where there is secondary emission nearly parallel to the initial partons. Several classes of ob- 
servables studied in Refs. [62-64] could be used to suppress exponentially forward radiation, 
which would then ensure the suppression of non-global, including super-leading, logarithms. 

In the following section, we study the definition and kinematics of rapidity gap processes 
and review the relevant factorization procedure. Here, we introduce the energy flow defined 
in Ref. [10], and define the gap fraction. After studying the gap fraction in terms of energy 
flow, we study its relation to color flow at short distances in Sees. [3] and HI We argue that the 
gap fraction exposes the pattern of soft gluon radiation, depending on the color SU(3) gauge 
content of the resonance in qq(gg) — > Z'/G — > QQ. In Sec. El the gap fraction for heavy 
quark pair production induced by a heavy spin-1 resonance is studied. Finally, in Sec. EJ we 
apply our formalism for a study of spin-0 or 2 resonance particles associated with different 
partonic processes. 

2 Heavy Quark Pair Rapidity Gap Cross Sections 
2.1 Definitions 

In this paper, we will focus on rapidity gap cross sections for heavy quark pair production 
through a new physics particle resonance denoted by V'(M,sp) of mass M and spin sp. 
The partonic scattering process for proton and (anti)-proton collisions, creating s-channel 
resonance V, is 



for the production of a heavy quark pair at fixed rapidity difference, Arj = i] a — rji,. Here 
the fi refer to partons that participate in the collision, and we have labeled process (pQ) "f" . 
We sum inclusively over final states, while measuring the energy flow, Q c , into the rapidity 
region, Q gap , in Eq. (CQ). The transverse momentum pt of a parton Q{p a ) or Q{pb) of mass 
rriQ in Eq. ([!]) is related to the rapidity difference At] according to 



f: /i(pi) +/ 2 (P2) - V'(M,sp) - Q(p a ) +Q( Pb ) +n gap (Q c )+X, 





(2) 
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As discussed in the introduction, we may expect different patterns of radiation from inter- 
mediate resonances if they are distinguished by different gauge content. 

The geometry for the gap at a collider is presented schematically in Fig. [TJ where the 
gap region is determined by rapidity range Y. The factorization that we will discuss in 
Sec. 12.21 enables us to work perturbatively. The partonic cross section at leading order for 
process (OQ) involves five external particles. The kinematics of this analysis are described in 
the Appendix. 



Figure 1: Light to heavy process with gluon radiation of momentum into a rapidity gap. 
The figure describes qq — > QQ process. Arj is the rapidity difference of a heavy quark pair, 
Q(p a ) and Q(pb), and Y is the rapidity range of the gap region. 1 and 2 denote the initial 
partons while a and b denote the final partons. 



2.2 Factorized and refactorized cross sections 

The inclusive cross section for heavy quark pair production at rapidity difference At] from a 
resonance at mass M with fixed energy flow Q c is given in standard factorized form by 

= £ Jdx 1 dx,^ /A (x uflF )^ /B {x i ,^) 7 ^ w (3) 
/1J2 

Here <\>f x iA and 4>/ 2 /b are parton distribution functions (PDFs), evaluated at the factorization 
scale fip. Below, we use the scale fip = px, the transverse momentum pt in Eq. (T2]). The 
partonic scattering cross section d&® /dAr)dQ c can be expanded in a s , starting from the 
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lowest order (LO) Born cross section 

da® da^ 



S(Q C ) + ■■■, (4) 



dAr) dQ c dAr] 

where corrections may include potentially large logarithms of Pt/Q c - We recall that the 
index f denotes the partonic process in Eq. ([1]). Following Ref. [11], we consider the partonic 
cross sections integrated over the transverse energy, Q c , radiated into a symmetric, central 
rapidity gap region, Y in Fig. [TJ up to a fixed value Q , 

da®(M,m Q ,Q ,]i F ,Ar],Y,a s ([i F ))_ f Qo ^ da® (M, m Q , Q c , ]i F , Arj, Y, a s (]i F )) ^ 



dAr] J Q dAr] dQ c 

When Eq. (EJ) is divided by the total cross section, the ratio can be interpreted as the 
probability for heavy quark pair production with soft gluon radiation into the gap region up 
to energy flow Q . For simplicity of notation, we keep dependence on the resonance width 
T implicit but recall that we treat Y/M as a number of order unity. 

The partonic cross sections of Eq. (jSJ) can be further refactorized into a hard scattering 



matrix H® and a soft matrix S®, 
da® 



dAr] 



(M,m Q ,Q ,]i F ,Ar],Y,a s (fi F )) = ^HP L (M,m Q , fi F , fi, Ar],a s {fi)) 



la 



xS®(Ari,Y,^,a s (fi),m Q j , (6) 

where matrix H)l represents interactions at short distances, independent of soft gluon radi- 
ation, and contains the dynamics of resonance V'(M, sp). Here we have introduced a new 
refactorization scale ji. The soft function S LI describes the radiation of soft gluons up to 
the scale Qq, which decouples from the dynamics of the hard scattering. The soft function 
is written in terms of products of path-ordered exponentials, w L (x)^}, of the gluon field 
as [11], 

S®= /^d^XXW^-tfJffl (7) 

n {h} 







where we sum over all states n whose transverse energy flow into the gap is restricted to 
equal Q' c and where color tensors labeled by L and / account for the color flow at the hard 
scattering. Eikonal multipoint operators, (wf^x))^}, with color indices are written as 
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products of ordered exponentials, tied together by color tensors, cl, 

{<k} 

x(cf) dMud2 <S> { £\0, -oo; x) dlifcl $^ 2) (0, -oo; x) d2 



&2 ) 



where the non-Abelian path-ordered phase operators (Wilson lines), $P , are given by 



&I\oo, 0; x) = Pexp 



-ig / d\(3-A^\\(3 + x) 







(9) 



Here, /3 is a four-velocity, and the vector potentials are in the color representation 
appropriate to flavor /. At the tree level, Eq. (jTj) reduces to the trace of the product of color 
tensors. 

The product of hard and soft functions in Eq. ([6]) is independent of the refactorization 
scale, p. From this, we easily derive for S^j the evolution equation [2] 



(10) 



in terms of a matrix of anomalous dimensions, (Tg)ji- The matrix (A77, Y,p) depends 
on kinematics, including the geometric information of the gap region and the heavy quark 
mass. Solving this equation will enable us to resum the leading logarithms of the soft scale 
Qq. In heavy quark pair production, dependence on quark mass and the partonic center of 
mass energy M in soft functions will always appear through the parameter p, given by 



^wi+r^V= , ! (id 



Pt 



1_ A cosh 2 4? 



1 M 2 2 

We will find that the soft function depends on the velocities (3i only through A77 and p. 



To obtain the solution to Eq. (flUj) . we treat above equation in a basis that diagonalizes 



r (f) 
1 s ' 



(rf(A V ,Y,p)) j/3 = \f(Ar),Y,p)8. 



/3 K^'h ^Pju-ifi 

= (R^UiTfiAvXp^uiR^U, (12) 
where the eigenvalues are given by a series in a s , 

\f(A v ,Y,p) = ^Xf 1 \A V ,Y,p) + -... (13) 

7T 
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We can transform the soft and hard matrices in the diagonalized basis 

[R%KH® L [R^) Lf3 . (14) 

( e /?)x) where ep are 



This change of basis uses the transformation matrix [(-R®) 1 }ki3 - 
eigenvectors of the soft anomalous dimension matrix. 

The solution to Eq. (fTDl) gives the dependence of on the ratio n/Q 



o- 



o(f) 



Ar/,^, — ,a a (fi),p 
P 



&fp(&ri, Y, 1, a.(Q ), p) exp 



Qo ^ 



2^ 



Sg(Atj,y,l,a.(Q ),p) 



7/3 



• (15) 



We will choose below // and \ip to be the hard scale pr for use in the refactorized partonic 



cross section 
of T 5 by 



In Eq. (j!5p the exponents i?^ are given in terms of sums of the eigenvalues 



E®(A v ,Y,p) = j^' 1 >(Arj,Y,p) + \f 1 \A V ,Y,p)} 



(16) 



where (3q is the lowest-order coefficient in the expansion of the QCD beta function, fi = 
(11N C — 2n/)/3, for N c colors and rif quark flavors. Equation (fl6l) shows that the diagonal 
exponents, are real. In the second form of Eq. ( fTBT) . we have used the one- loop QCD 
running coupling, 

2tt 



(17) 



A)ln(/i/A)" 

Combining Eqs. and (|15|) . the partonic cross section, valid to leading logarithm, is then 



dAr] 



J2H^ LO \M,m Q ,Ar),a s ( PT ))S 7 p 



(f,o) 



0,7 



7/3 



i / Pt 



? (f) 



18) 



In the following section, we will use Eq. ( fl8i) to evaluate the gap fraction. 



2.3 The gap fraction 

We define the gap fraction f gap as the ratio of the number of events for heavy quark pair 
production with a specified rapidity gap to the total number of the pair production events. 
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A gap event was originally identified experimentally by the lack of particle multiplicity in 
the interjet region [5,8]. In that case, the multiplicity is determined from the number of 
calorimeter cells which measure energy deposition above a threshold. A gap event is then 
defined by the absence of such cells in the rapidity region. Our formulation of the problem is 
in terms of the transverse energy flow, Q c , of hadronic radiation carried by the particles [10]. 
We introduce a variable, energy "gap threshold" Qq, which is different in principle from the 
experimental calorimeter threshold, and identify a gap event from the condition of interjet 
radiation less than Qq. Then, to get the partonic gap cross section, we evaluate Eq. (J3J) at 
the value Qq for the maximum of interjet energy flow Q c . The resummed result is given 
by Eq. (TT8l . We note that values of Qq at the order of the hard scale pt would violate the 
requirement of factorization, since the emission into the gap region would not be soft any 
more. 

To define the gap fraction, we approximate the total cross sections from the LO (Born) 
partonic cross section, related to Eq. (jSJ) 

^ = Y, H lt°\ M ^^(p T ))S^. (19) 

We may consider the LO and the gap differential partonic cross sections with respect to M 
and r)y = | In given for a resonance by, respectively, 

rfgftso) ( M , m q> A V , a a (j>r),Y) do^ ( 2 . ( 1 x x 

— — f — — = — — o [M — XiXob ) o rpji In — 

dAr}dM 2 dr} V , dAr] v ; V lv 2 x 2 

do® (M, m Q , Q , Atj, a s (p T ), Y) da® . ( 1 xA 

dA v dM*d Vv , = W ( 1X2 ] 5 { VV '~2 ln x- 2 ) 



(20) 



which are distinguished by the argument Qq in the latter, and where S is the hadronic center 
of mass energy. Here, the differential with respect to rjy corresponds to fixing the center of 
mass rapidity of the final-state pair in the lab frame. We emphasize that the above condition 
does not change the description of soft QCD radiation in the soft functions, which depend 
on the rapidity difference Ar). 

We now define the gap fraction at leading order in terms of Eqs. ( TTBl . (|T9|) and (1201) . by 



(LO) = 22 hJ J dx 1 dx 2 (f )fl/A {x l ) ( j )h/B (x 2 ) d ^ dM -d Vv , 

hap t »(^O)(M ) m ,A^a i (pT),y) ' [ ' 

E /ll/2 J dx x dx 2 <t> hlA {x x )<j> hfB {x 2 ) y d A V dAPd Vv , ' 
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An advantage of using the differential partonic cross sections (|20|) to obtain the gap fraction 
becomes clear when there is only one partonic initial state for a resonance process. The gap 
fraction at leading order in Eq. ( l2TT) for a single partonic channel process is given by 

doM (M, m Q , Qq, Ar,, a a (p T ), Y) 

&> W m Q , q„, A,, a.(Rr), K, = Ja ,»o, (M , m f !„,„,, „,, y y . (22) 
where the PDFs cancel in the ratios of the gap fraction in this case. 



3 Massive Gap Soft Anomalous Dimension Matrices 

The one-loop soft anomalous dimension matrix, Ts, for a rapidity gap process is calculated 
from the mismatch between real and virtual corrections, generated by imposing a rapidity 
gap. The soft function is calculated by taking single ultraviolet pole parts from virtual and 
real corrections, which we will label as u)\r(ij) an d UJ R(ij)i respectively. The real part, ujru^, 
comes with a phase space integral outside the gap, where it cancels uJvtij)- The coefficient of 
the pole results from the phase space integration inside the rapidity gap region. Here, the 
single ultraviolet pole arises from the limitation on the energy integration of soft gluons. 

In practice, following [11], we construct the gap soft anomalous dimension matrix from 
a set of integrals, u>(ij), 

(1) _ ^ fa -)(-2e) 

where we use dimensional regularization, with D = 4 — 2s. The indices i and j label partons 
of momenta Pi and pj, and the 0J{ij) are defined by 

-(47ra.)<WA4 1 ^ x iPl ' Pj) 



+5 i 5 j A i A j ^(l-5 i S j ). (24) 

Here the subscript P.P. indicates that the integral is defined by its ultraviolet pole part. 
We also define 5, = 1(— 1) for momentum k flowing in the same (opposite) direction as the 
momentum flow of line i, and Aj = 1(— 1) for i a quark (antiquark) line. The function Q(k) = 
1 when the vector k is directed into a rapidity gap. Equation (|24"|) describes the incomplete 
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cancellation between real and virtual gluons from the integration over the geometric phase 
space of the gap region. Only the part of the phase integration inside the gap survives, and 
contributes to the soft anomalous dimension matrix for the process. 

We express the above integrals in terms of transverse momentum, rapidity, and azimuthal 
angle coordinates, with the help of the kinematics described in the Appendix, and use the 
relations 



d D k5 + (k 2 



p.p. 



h.D-3 



dkxdydcj) . 



p.p. 



h 1+2s 

P.P ""7 1 



1 

2e 



(25) 



In these terms, the u^j) become 



1 1 



3 3 2 2e / 2tt y ' 3 3 3 2e 2 



(26) 



where we define fly as 



(Pi ■ k)(pj ■ k) 
In these coordinates, the fl^- are given by 

Sli2 = 2 , 



(27) 



p 2 cosh Arj + 1 



ab 



n 



ib 



n 



2d 



2b 



(pcosh (^ — y) — cosh0 


) (p cosh (4r + y) + cos <p) 


pe-(^-y) 




(pcosh {pjr — y) — cos0) 


) 


pe(^) 




(pcosh {^2~ + y) + COS0) 


? 


pe^- y ) 




(pcosh (^ — y) — cos0) 


i 






(p cosh (^ + y) + cos (fj 





(28) 



We recall that indices 1 and 2 are for initial-state partons, and a and b for final-state partons, 
as described in Eq. (JT]). The parameter p has been defined in Eq. f ill I) , and is always larger 
than one for massive partons and equal to one for the massless case. We will see that 
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all the results we derive here go back to the previous calculation for massless dijet events, 
p=l [10,12]. 



For complete integrations over the rapidity gap geometry, it is convenient to integrate 
Qij over an azimuthal angle <fr first. It is then straightforward to evaluate rapidity integrals 



over a gap region — Y/2 to Y/2. The final expressions for the Isuj-j of Eq. ( J23i) are then 



J?) 
IS (12) 

J 1 ) 
IS (ab) 



x In 



IS (la) 



Y — iir 
—in + 



p 2 cosh Ar] + 1 



2p cosh(^) Jl + p 2 sinh 2 



cosh (A^) sinh (*z) + sinh (Z) + J( p2 sinh 2 (^z) + X ) ( p2 cosh 2 (A^yj _ ^ 



p 2 cosh (AztT) sinh ( A2) _ sinh g) + J( p 2 sinh 2 (A,) + X ) ( p2 cosh 2 (A^) _ ^ 



111 



(P 2 + 1) 



2p (pcosh 2 (^) + sinh (^±L) ^Jp 2 cosh 2 - 

- (p 2 + 1) 



/ 



y 2p (^cosh 2 (^) + sinh yV cosh 2 (^1) 



7s (16) 



-iln 
2 



-iln 

4 



pcosh (^) + Jp 2 cosh 2 - 1 



pcosh ( Av 2 Y ) + yp 2 cosh 2 ( A?? 2 y ) — 1 
^2p (pcosh 2 + sinh ^h^f^ij - (p 2 + 1)^ 

^2p (pcosh 2 (*t£) + sinh (£t£) ^p^rfT^Pl) - (P 2 + 1) 



pcosh (Y) + ^sh 2 (^)~, 



'S (2a) 
IS (2b) 



— In - 

2 I 

IS (lb) • 

IS (la) 



(29) 



We take the massless limit, p = 1, in Eq. (I2"9~|) to cross-check these results with the previous 
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calculation for massless dijet events, and find [12] 

Y — m , 



7S (12) 
IS (ab) 



hi 



' sinh (*H±L) 



sinh(^r) 

!? /sinh(^) 

75 (la) " 75 (2b ) " 2 (^ ^^(A^H) ' " V 



7 « - 7 f) lA/ sinh^) 

We can now construct the soft anomalous dimension matrices which depend on the color 
basis. We use an s-channel singlet-octet basis [51] for quark- ant iquark annihilation to a 
heavy quark pair, 

Cginglet = c li Coctet = Cl ^ ' (3-Q 

where c\ and C2 are color tensor basis for singlet exchange in the s and t channels, 

( C l){n} = firir 2 3r a r b y {C2){ri} = ^nra^nn ■ (32) 

Here, we denote by fj the color index associated with the incoming or outgoing parton i. The 
lowest-order soft matrix is the trace of the color basis, S^f* = Tr (c^c/), given in singlet-octet 
basis by 

/ N 2 C \ 

= , (33) 

V \{N- c -l) ) 

which will be used in the resummed cross sections, Eq. ( IT8|) . 

To construct the soft anomalous dimension matrix in the color basis (j3ip . we rewrite the 
various one-loop diagrams, using the identity shown in Fig. |2j 

( T F)ji( T F)ki = ^ (tikiSji - jj-SjiSiaJ , (34) 

where the T^'s are the generators of SU(iV c ) in the fundamental representation. In the basis 
of Eq. ( 13~TT) for the process qq — > QQ, we derive in this way the soft anomalous dimension 
matrix in terms of the integrals l^dj) m Eq- ( |29l) . 

r a )( ^QQ) = ( ^ - ■ - j (3rj) 




Figure 2: Color identity corresponding to Eq. ((3 
where a, X an d r are defined by 



a 


— IS (12) + 7S (06) ' 




X 


— '5 (la) + 75 (26) ' 




T 


— 75 (16) + 75 (2a) ' 


(36) 



The eigenvectors of in Eq. fl35l) may be chosen as 

( = j 47V c g 



1 

and the corresponding eigenvalues are given by 



e 2 = I InTz I , (37) 



(fil) A 2 + ^A? + 8C F JV C £ 2 
A i - 



4iV c 



■ 



A (f J i) = A 2 - ygj + 8C F N c e m 

where we have defined 

Ax = {a + x + 2T) + 2C F N c (a- X ), 
A 2 ee -(a + x + 2r) + 2C F N c (a + X ) , 
£ = (r + X ). (39) 

Finally, the transformation matrix R in Eq. (1121) . which diagonalizes the soft anomalous 
dimension matrix, is given from Eq. (13"T|) by 



(^ (fp1 )^ = ( *L A Y I • ( 4 o) 

Ai+t/A^+SOA^ 2 
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It was observed in Ref. [10] that for massless partons the above eigenvectors are independent 
of the rapidity difference Arj, and depend only on the rapidity gap range, Y. However, we 
find that this is not the case any more for massive quark pair production, p > 1. As an 
example, we take the values Arj = 2.5 and At] = 4 at fixed mq = m t , M = 1.5 TeV and 
Y — 1.5, and confirm a mild rapidity dependence of the eigenvectors, 

e ' = ( -0.098 - 0.269i ) forA 'f = 2 - 5 . e ' =(-0.102 - 0.274i) forA '' = 4 ' 
e 2 = ( °' 022 + 060i ) for A, = 2.5, e 2 = ( °' 023 + a08M ) for A, = 4 . (41) 

For these configurations, the first eigenvector, e±, is close to a color singlet, while the second, 
e2, is close to a color octet. Following the terminology in Ref. [10], we denote the eigenvectors, 
e\ and e 2 , of T s by "quasi-singlet" and "quasi-octet" , respectively. 

We will see below that the process qq — > QQ is the only relevant partonic channel 
involving a spin-1 particle resonance, V'(M,sp = 1). For other resonances V'(M,sp 7^ 1), 
we can consider the gluon-induced partonic process 

9(pi) + 9(P2) V'(M, sp 1) ^ Q(p ) + Q(p b ) + % ap (g c ) . (42) 

To describe the color flow in the latter process, a suitable s-channel basis of color tensors 
has been defined in Ref. [11], 

{ c 2){n} = d rir2C (Tp) raTh , 

( C 3){rJ = ifr ir 2c(Tp) rar . b , (43) 

where c\ is again the s-channel singlet tensor, and C2 and C3 are the symmetric and anti- 
symmetric octet tensors, respectively. In this basis, the anomalous dimension matrix for the 
process gg — > QQ is found to be [12] 

f C F ^\ ab) +N c 4\ l2) |( X + r) \ 

iV c C f(x + r) , (44) 

V X + r ^nt(x + t) N c C J 



r (i)(9S-QQ) 
1 5 



where £ = |(% — r) — ^^Ts ( a &) + I^S (12)' anc ^ where x and r have been defined in Eq. f[3T)j) . 
In analogy to Eq. (|33l) . we derive the zeroth order soft matrix for this process in basis (1431) 
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as 



5-Cf.o) 





V o 



(jVg-4)(jVg-l) 








2 



(45) 



/ 



One may now find the corresponding eigenvalues and transformation matrices 1 for 
99 ~~ * QQ- The results are too cumbersome to present here, but are straightforward to derive 
and evaluate numerically. In large N c limit, however, the eigenvectors and the corresponding 
eigenvalues can be obtained readily, since Eq. (J33J) becomes 



2 IS (ab) + IS (12) 



V 



for large N c . Normalized eigenvectors of Tg 



ei 




e 2 




(46) 



(47) 



and the corresponding eigenvalues are given by 



(f>i) 



(i) 



7sV) +7s'(i2) J 



A 



(f>i) 



C--( X + r))JV c . 



(48) 



We emphasize that, in Eq. (j47p . ei is an exact color singlet and that e<i and e3 are octet 
basis vectors. We will not restrict ourselves the above approximation in following sections. 
We will see, however, in Sec. |6] that the explicit results of eigenvalues and eigenvectors for 
99 ~ > QQ i n the kinematics and the gap geometry we study in this paper are close to the 
results of the large N c limit. 

The analysis in this section is based on wide resonances with r ~ M of decay width, 
T, and resonance mass, M. We note, however, that the energy scale fi of gluon radiation 
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into a gap should be separated into two different regions, depending on the decay widths of 
heavy resonances. This may lead to different prescriptions for describing soft gluon radiation, 
depending on fi. For gluon radiation at scale //, with Q$ < ji < VTM, the hard function can 
be treated as an effective vertex. This is the condition we have used for describing soft gluon 
radiation in this section, where we assume that the hard scale px is of the same order as 
VTM. For the case of gluon radiation in the range < H < M, direct radiation from the 

resonance, which was not considered in this section, has to be included in the soft function. 
In addition to this, the contributions from soft gluon interactions between initial and final 
state partons are suppressed for both octet and singlet resonance processes if the hard scale 
falls into the regime \/YM < jj < M. In this paper, we concentrate on the case that VTM is 
of the order of PtS I n general, the decay width of an octet resonance would be large enough 
to satisfy this condition. For example, a KK gluon above 1 TeV (as required by precision 
tests) has decay width of about Mg/6 [27,29] in basic RS models, and 0.215Mg [32,65] in 
RS models based on extended electroweak gauge symmetry with specific fermion charges 
and localization that explain the LEP anomaly of the forward-backwardy asymmetry for 
b quarks. Axigluons or universal colorons have decay width of about a s Ma [48], resulting 
in y/TMc ~ Pt where pr is the hard scale we have used. On the other hand, the widths 
of singlet resonances are sensitive to the model parameters for the Z' . A Z' can be very 
narrow, O(10~ 3 Mz'), or very broad, O(Mz'), in both Topcolor Z' models [66] and in the 
Little Higgs models [67-70], depending on parameters such as the S77(2)i <g> SU{2) 2 mixing 
angle. In the following sections, our study will be based primarily on the broad widths of 
heavy resonances, Qq < fi < yVM. We have observed that for narrow resonances we have 
to evolve /i both from Qq to yTM with the exponents we have obtained in this section and 
from \/YM to the hard scale with a different analysis. We will discuss briefly how such a 
narrow singlet resonance can change our results at the end of Sec. |6j 

4 Hard Functions 

The remaining piece to complete constructing a partonic cross section in Eq. ([6]) is a hard 
function. Hard functions describe physics at the large momentum scale. In our study, 
the mechanism of new heavy particle resonances is thus contained in the hard functions, 



1 In principle, a hard scale for the specific model could be chosen as \/TM , the minimum off-shellness of 
the resonance. The hard scale pr is suitable in our study since we are interested in a model independent 
approach and since we consider broad resonances that satisfy ^/TM ~ px by assumption. 
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decoupled from soft QCD radiation. The direct radiation from an octet resonance can thus 
be neglected. In the previous section, we have constructed soft functions in the space of 
color flow at short distances. Thus, it is convenient to construct hard functions in the same 
basis. In Fig. [3], we represent several possible diagrams for the hard part, involving s-channel 
heavy resonances at lowest order. In the color basis (I3T1) . the corresponding hard functions 
for qq — > G / Z' — > QQ in Figs. 3(a) and 3(b) are given by, respectively, 



[H, 



{t,Loy 

G 



IKL 



h G (M,s,t,u) 




1 



(49) 



IKL 



hz'(M, s, t, u) 



1 




(50) 



A hard function for gg 
would take the form 



G — ► QQ in Fig. 3(c) such as a KK gluon resonance process [44,45] 





(H'^ LO) ) KL = h' G (M,s,t,u) [ 

1 



(51) 



and for gg — > Z' — > QQ in Fig. 3(d) a hard function is given by 



'( f . L °)\ 



IKL 



h' z ,(M, s, t, u 




(52) 



in the color basis (j4"3"j) . The overall functions ha/z' (M, s, t, u) and h! G , z , (M, s, t, u) depend on 
hard scale kinematic variables, M, s, t and u, independent of the gap variables, Qq and Y . 
We note that the spin of the particle V limits possible partonic cross sections. Specifically, 
as noted in Refs. [71,72], the Landau- Yang theorem [73] forbids the creation of a massive 
state of total angular momentum J = 1 from two massless spin-1 particles. Therefore, the 
gluon-induced production of a spin-1 massive state in Fig. 3(c) or 3(d) is not allowed if the 



spin of the resonance is one. Spin-0 singlet or octet resonances could arise, however, in 
split SUSY scenarios via gg annihilation. In such cases, gg — ► G/Z'(sp = 0) could be the 
dominant production channel, giving hard functions of the forms, Eqs. (j5ip and fl52|) . Also, 
spin-2 particles V'(sp = 2) such as gravitons and reggeons [74] could be created in s-channel 
resonances through gg — > V'(sp — 2) — > QQ, as well as qq — > V'(sp — 2) — > QQ. The cross 
sections for resonances in these cases depend crucially on the model we study. For spin-1 
resonance processes, we do not need to know specific coupling strength to the SM particles, 
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because spin-1 G and Z' involve only the quark- induced partonic process, qq — > G/Z' — > QQ. 
This simplifies the study of the gap fraction. According to the argument in Sec. I2.3[ the gap 
fraction at leading logarithm results in the ratio of partonic cross sections ( 122|) . free from 
the PDFs that cancel in the ratios. Also, ha or hz of the hard-scattering functions cancels 
in the ratios ff22|) . We can thus study the gap fraction analytically. 

In Sec. El we will discuss how gluon-induced resonances contribute to the gap fraction in 
massive spin-0 or 2 resonance processes. For a complete study of the gap fraction in this case, 
we should combine all the possible partonic processes that depend on the specific model. In 
addition to these, we must convolute the partonic level results with PDFs, following Eq. ([3]). 
We leave this study for a further project. 

To obtain the partonic cross section in Eq. ffl8l) . hard functions in Eqs. fj49l) . fl50l) . fl5Tj) . 
and (1521) are transformed to (H^' LOS) )^p through Eq. (Tl4j) in the color basis that diagonalizes 
T s in Eq. ( TT2l) . We will study explicitly how partonic cross sections and gap fractions are 
evaluated in the following sections. 
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5 Gap Fractions for Spin-1 Resonances 



From the previous sections, we have all the tools necessary for the evaluation of the partonic 
cross section in Eq. ( |T8l) . As mentioned above, the main quantity we focus on is the gap 
fraction (J2TT) . For spin-1 resonances, we can limit ourselves to quark- induced resonances, 



shown in Figs. 3(a) and 3(b) , following the argument in Sec. HI The partonic cross section 



for heavy quark production via G or Z' resonances is 



da 



(f) 

G/Z' 



E(4 f }z 0) )/37(M, mo, A77, a s (p T )) S« 0) 



dAr] 



J 7/3 



(53) 



where the (H^^)^ are hard functions written in the basis that diagonalizes the soft anoma- 



lous dimension matrix. We will implement Eq. fl53l in numerical studies below. Following 
Eq. (fill) , to obtain (Hq^)^, we transform the hard functions in Eqs. ( pf9l) and ( |50l) . 
through matrix R given explicitly in Eq. ( )40l) . The soft matrix S^g is obtained by trans- 



forming Eq. (133]) through the same matrix R. Finally, the exponents El that appear in the 
cross sections (153"]) are given by Eq. (|T6|) . 



,(f) 



In Eq. f l53|) . we observe that the Qo-dependent factor, [in i^f) / In (^f )] 7/3 , is larger at 
fixed Qo for smaller E^, because [in (^) / In (^f)] is less than one for Q < pr- Note that 
we apply Eq. ( 153]) only for Qo < Pt- A constraint for the maximum value of Arj is obtained 
from p in Eq. ( TTTT) . 



y < A77 < 2cosh _1 ( ) . 



(54) 



Recall that the gap fraction measures the fraction of events for which the radiated trans- 
verse energy is less than or equal to the gap threshold Q . For any E^ > 0, the fraction 
vanishes for Q = A, and approaches unity for Q = Pt, where we choose /i = pt in numeri- 
cal calculations. The soft anomalous dimension matrix for this analysis is given by Eq. (J35J . 
Representative numerical values of the real and imaginary parts of the eigenvalues are shown 
in Fig. H] as a function of gap size Y . The corresponding exponents at fixed Y are shown in 
Fig. [5j as a function of At]. Only the real parts enter the cross section. Imaginary parts are 
shown for completeness. Qualitatively, in Figs. @] and El we see that the eigenvalues and the 
corresponding exponents are larger for larger Y and for smaller At]. In Fig. HI we observe 
that for y of order unity, Re[Ax] is much larger than Re[A 2 ], leading to E$ < 1 < E§, as 
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(a) (b) 



Figure 4: Plot of the real (a) and imaginary (b) parts of the eigenvalues of the soft anomalous 
dimension matrix in Eq. (j35p for qq — > QQ, as a function of Y for M = 1.5 TeV, At] = 2.5 
and rriQ = m t . The solid line identifies the quasi-singlet eigenvalue, A^ , the dashed line, 
quasi-octet, A| , in Eq. (137|) . 




Figure 5: Plot of the exponents, -Eia, of the soft anomalous dimension matrix for qq — > QQ, 
as a function of Arj for Y = 1 and mg = m t with resonance mass M = 1.5 TeV and M = 
750 GeV for (a) and (b), respectively. The solid lines identify the quasi-singlet exponent, 
E®, the dashed lines, quasi-octet, E®, obtained from Eq. ffT6]) . 
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shown in Fig. [5j Following the discussion below Eq. (|53|) . a smaller value of the exponent 

i?22 enhances the factor [in ( Q ) / In 22 m the partonic cross section (I53I) relative to a 

larger value. We note that off-diagnal contributions are small. 

The resulting gap fractions for two different spin-1 resonances, Z' and G, of mass 2 TeV, 
1.5 TeV, 750 GeV, and 550 GeV are illustrated as functions of Qo in Fig. [6] as the solid and 
the dashed curves for a top pair, and as the dot-dashed and the dotted curves (blue) for a 
bottom pair at fixed Y — 1.5 and Arj = 2. The ratios of gap fractions, f^ ap /fga P , are shown 
in Fig. [7] at fixed Q = 5 GeV in terms of Y and Ar] for M = 1.5 TeV and 750 GeV. We 
observe that for the ratios the dependence on At] is negligible at fixed Y. Also, Fig. [7] shows 
that these ratios get larger at fixed Qo as the gap range increases. 

In Figs. M and [HI the gap fractions and their ratios are presented as functions of Q and 
Y for M — 1.5 TeV and 750 GeV, respectively. The gap fraction for an octet resonance 
increases rapidly for small values of Qo i n comparison to a singlet resonance for any gap size. 
Many fewer events are accumulated at low Qo through a color-singlet resonance than for a 



color-octet resonance. For instance, at Qo = 3 GeV and V'(1.5 TeV, sp = 1) in Fig. 6(b) 
we obtain gap fractions of 12% and 54% for color-singlet and -octet resonances, respectively. 
The figures also show that the ratios are larger for larger M and Y. 

We can interpret these results as follows. As we saw in Eq. ( ]411 of Sec. El the eigenvectors, 
ei and e2, of the soft anomalous dimension matrix for qq — > QQ are very close to a color 
singlet and a color octet for the kinematics and the gap geometry we study in this section. For 
this reason, we have called them "quasi-singlet" and "quasi-octet" eigenvectors, respectively. 
Consequently, this leads the hard functions in Eqs. ( 1491 and ( 1501 to retain approximately 
their color basis dependence even after transforming the functions to the diagonal basis. 



Following the above property of the hard functions and the eigenvalues, shown in Fig. 4 (a 



for qq — > QQ, the quasi-octet component of the hard function for a color octet resonance is 

(f) 

the leading component, with a corresponding small exponent, E 22 . We thus find that the 

Q E 22 

quasi-octet cross section (H^' L °^) 2 2S22°^ rjwl is dominant in the partonic gap cross 
section (153!) for an octet resonance G. The quasi-octet cross section is enhanced due to its 
large overlap with the color-octet basis which gives a small exponent E22 < 1. It is clear 
that this term produces the "convex" shapes of the dashed gap fraction curves (f g ap) of the 
log plot on the shown in Fig. [6J and the upper surfaces in Figs. [8] and M for top 
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M= 2 TeV 



M= 1.5 TeV 




Q 



<» 



lGeV 




(a) 

M= 750 GeV 



(b) 

M= 550 GeV 







lGeV 




50.0 100.0 j Q e y 



Figure 6: The fractions for gaps identified by the energy threshold Qo at Ar/ = 2 and 
Y = 1.5 for resonance masses M — 2, 1.5, 0.75, and 0.55 TeV. In the above figures, the solid 
curves describe the gap fraction through a Z' resonance (color- singlet), and the dashed curves 
throughout a G resonance (color-octet) into a top quark pair. The dot-dashed curves (blue) 
describe the gap fraction through a Z' resonance (color- singlet), and the dotted curves (blue) 
throughout a G resonance (color-octet) into a bottom pair. 
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(a) (b) 

Figure 7: The ratios of two different gap fractions, J fg , as functions of Y and Ar/ at 
fixed Qq = 5 GeV are shown (a) for M — 1.5 TeV and (b) M = 750 GeV, respectively. 



pair production. These curves rise rapidly at low Q 



o- 



The singlet resonance process, on the other hand, has a larger overlap with the quasi- 
singlet component of the hard function. The eigenvalue of the quasi-singlet is larger than one, 

m(%) 



and -E^i = (4//3o)ReAi > 1. Thus, the quasi-singlet cross section, (H^', L0 ^)iiSil 
is suppressed at low Q , explaining the "concave" shapes of the solid lines in Fig. [of and the 
lower surfaces in Figs. [Hand Of or top pair production. These increase slowly at low Qq. The 
differential partonic cross sections with respect to Q c are presented in Fig. [10] for octet and 
singlet resonance processes for values M = 1.5 TeV, Arj = 2.5, and Y = 1.5. The behavior 
of the exponents corresponds to the dominance of the octet resonance process at low Q c in 
Fig.EOj 

The gap fraction depends on the quark mass, ttlq, through the variable p in Eq. (TTTj) . 
The dependence can be seen in Fig. O by comparing the black curves (solid and dashed) for 
a top pair with the blue curves (dotted and dot-dashed) for a bottom pair. For the b quark, 
the gap fractions are slightly smaller at fixed gap threshold energy than for a top quark pair 
for the relatively heavy resonances in Fig. [6j This difference is, however, larger for smaller 
resonance mass, as also shown in Fig. [HI We note that the quark mass dependence on the gap 
fraction for bottom quark pair production is negligible even for a resonance of relatively low 
mass, but not for a top pair. For instance, we obtain p = 1.0003 for bottom pair production 
via a resonance of M — 550 GeV with Arj = 2, while p = 3.8200 for a top pair. This results 



in the very different shapes of the black curves compared to the blue curves in Fig. 6(d) 
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(c) (d) 

Figure 8: Gap fractions as functions of energy threshold Qq and gap range Y from a 
resonance of M — 1.5 TeV, (a) for At] = 3.5 and (b) for Arj = 2.5. The lower surfaces in (a) 
and (b) describe the gap fractions through a Z' resonance (color-singlet), the upper surfaces 
throughout a G resonance (color-octet) decaying into a top quark pair. The ratios of the 
gap fraction for an octet resonance to the gap fraction for a singlet resonance are illustrated 
in (c) and (d) as functions of Q and Y. 
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Figure 9: Same as Fig. [HI for M = 750 GeV, in the allowed region in gap size Y. 
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Figure 10: The cross section from the octet resonance (dashed line) and the singlet res- 
onance (solid line) for M = 1.5 TeV, At] = 2.5, and Y = 1.5. Both distributions are 
normalized to the same area for the range of A < Q c < px- 

The results in this section show that the measurement of gap fractions can be an effective 
tool to identify the color content of a spin-1 s-channel resonance. In the next section, we 
study gap fractions for spin-0 or 2 resonance processes, which contain more partonic channels. 



6 Gap Fractions for Spin-0 or 2 Resonances 

As we discussed briefly in Sec. HI in several models of new physics, there are s-channel spin-0 
and spin-2 resonances, decaying into a heavy quark pair, which can be produced, through 
gluon-induced processes in addition to quark-antiquark reactions. We readily recognize that 
spin-0 or 2 resonance processes might lead to a different pattern of soft gluon radiation 
between the Tevatron and the LHC, because of the presence of both quark and gluon- 
induced partonic processes. As we shall see, PDFs no longer cancel in the leading order gap 
fraction. This is distinguished from the analysis in Sec. [5] for spin-1 resonances, involving 
only one partonic channel. The results in Sec. [5] are energy- independent, at leading order, 
applying to both the Tevatron and the LHC. 

Let's assume that there are two partonic processes that produce a heavy resonance, 
decaying into a heavy quark pair, : qq — > V — > QQ and f 2 : gg —* V — > QQ. In addition, 
we assume for simplicity that the coupling of light quarks is universal. Recall that the process 
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f 2 cannot occur for a spin-1 resonance in general. Therefore, we consider the octet resonance 
as spin-0 or 2. In this case, we obtain the gap fraction, following Eq. (12ip . 

f(LO) 
J gap 



Jd Xl dx 2 (J2f lf2 dSfSt^, Hi/ A (Xl) <t>h/B {X 2 ) + d £f d % Q ^ v , <\>9lA {X X ) (pg/B (x 2 ) 



Jd Xl dx 2 (j2 fih 1 0^- ( f >fl/A (x 1 ) ( p h/B {x 2 ) + J 0^- ( p g/A ( Xl )<P g/B {x 2 ] 
pqq^QQ jqq^QQ _|_ pgg~*QQ jgg^QQ ^55^ 



i gap ' f J gap 

where the p fi are defined by 

p qq-+QQ = dA H 



Y. hh ( Ph/A{Xl)<P h /B{x 2 )^ 



v99 ^QQ = ggW^g^ m 



dAv (Ef lh <t>h/A{X\)4>h/B{x2) ) + a d l v <Pf g /A{xi)<j)f g /B{x 2 ) 

and the gap fractions, fgl^ QQ and fga^ QQ , are given by 

Jgap WO) - dtVi.LO) ' J 9ap WOj ~ d -(f 2 ,LO) > \° ' ) 

dArj dAr) 

which can be calculated perturbatively, following Eqs. ( 1T81) and (|T9l) . We can interpret the 
p fj as the probability for the process fj among the channels, p 9 ^<9<9 _|_ ^gg-^QQ _ \ 

We can study the above gap fraction for representative choices of p99-»QQ anc i pgg^QQ _ p or 
the first example, we compare the gap fraction for an octet resonance with the gap fraction 
for a singlet resonance through a pure gluon-induced process, p99^QQ = \ anc [ pqq-*QQ = o. A 
pseudo-scalar or a boson-phobic scalar resonance is a good example for the above resonance 
process, where gluon-induced processes are dominant and the branching ratio to a top pair 
can be taken to be unity [26]. In general, SUSY models with this feature can be constructed 
for singlet resonances [75], and split SUSY scenarios can explain octet resonances from bound 
states of meta-stable gluinos [76]. 

The gap fractions for the pure gluon-induced resonances are determined by the soft 
anomalous dimension matrix in Eq. ()44l) . its eigenvalues and eigenvectors for the process 
99 ~ > QQ- The resulting gap fractions for octet and singlet resonances are illustrated in 



Fig. 11(a) for M = 1.5 TeV, Arj = 2.5, and Y = 1. 



To study the effect of an admixture of quark- and gluon-induced processes, the gap 
fraction for an octet resonance with pdg^QQ — q.8 and p^9-»QQ = q 2 is represented in 
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Fig. 11(b) as a dashed curve in comparison to the gap fraction for the pure gluon- induced 
singlet resonance, shown as a solid curve in the same figure. This quark/gluon admixture 
for an octet resonance can occur for the reggeon resonances of Ref. [74] . 



M= 1.5 TeV 



M= 1.5 TeV 
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lGeV 




go 
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Figure 11: Gap fraction as a function of gap energy threshold Qo at At] = 2.5 with M = 
1.5 TeV and Y = 1. In (a) and (b), the solid lines describe the gap fraction through a 
Z' resonance (color-singlet) for p99^QQ — \ The dashed line describes the gap fraction 
through a G resonance (color-octet) in (a) for = 1 and in (b) for p99-*QQ = o.8 and 

pgq^QQ = o.2. 
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that there is "less" radiation 



As observed in the previous section, we again see in Fig. 
into the gap region for the color-octet resonance process. At the same time, we notice that 
for the same resonances there is "more" radiation for gluon-induced processes than for quark- 
induced processes. This can be seen by comparing solid and dashed curves in Fig. [6] with 



the curves in Fig. 11 (a 



Scattered gluons tend to emit more radiation, compared to scattered quarks. We note 



that the shape of the gap fraction for gluon-induced octet resonances in Fig. 11(a) is "linear", 
while the shape for quark-induced octet resonances in Fig. [6] is "convex" . Following the 
arguments in the previous section, the values of the exponents for gluon-induced processes 
in Fig. [12] can account for the difference. The eigenvalues for the process gg — > QQ in Fig. 



ensure that all exponents E^ 99 ~*^^ are larger than one for Y of order unity, as shown in 
Fig. d2J We note that two quasi-octet eigenvectors, e 2 and e^, of in Eq. (jHJ), which 

we recall are determined numerically, turn out to approach e 2 and e 3 in Eq. (JJ7J), obatined 
in large N c limit. As a result, both the quasi-octet components, (H' G )22 and (H'q)^, of 
the octet hard function, written in the basis that diagonalizes rW( fl, fl- >< 3 < 9)^ arg i m p 0r t an t. 
Their exponents, E^ 9 ^^ and E^ 9 3 9 ^^\ are shown as the dashed and the dot-dashed lines 
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1 1 ) — 

Figure 12: Plot of the exponents, Eaa, of the soft anomalous dimension matrix for gg — > QQ 
as a function of A?] for Y — 1.5 and m = ra t with resonance mass M = 1.5 TeV. The solid 
line identifies the quasi-singlet exponent, the dashed and the dot-dashed lines the two quasi- 
octets. 




Figure 13: Plot of the real (a) and imaginary (b) parts of the eigenvalues of the soft 
anomalous dimension matrix for gg — > QQ of M — 1.5 TeV, A?7 = 2.5 and tuq = m t . The 
solid line identifies the quasi-singlet eigenvalue, the dashed and the dot-dashed lines the two 
quasi-octets. 
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in Fig. [12j An interesting feature is that one of the quasi-octet exponents is close to the 
exponent for a quasi-singlet. A semi- numerical expression for the cross section (fl8l) for 
gluon- induced octet resonances is given for M = 1.5 TeV, Y = 1, and Ar) = 1.5 by 



da, 



(gg^QQ) 
G 



E m '(f 2 ,LO)x q 



dAr] 



(fc.o) 



0,7 



7 W 



h' G (0.23) 



MS) 
Wf) 



Mf) 

1 2.27 

+ (6.19) 



(5? 



-i 2.08 



+ (6.33) 



Ml) 



1.17> 



where the first coefficient, 0.23, of second line in Eq. 

T '(h,LO)s _ c (f 2 ,0) , (rr'{h,LO) 



is (^ f2 ' iO) ) n ^ 2 ' 0) , and the 



other two, 6.19 and 6.33, are (H G I2 '^ U) ) 2 2S22 ' u> and (^? I2, ^033^33 ' respectively. Here 



h! G has been defined in Eq. (l5Tj) . In Eq. (1581) . we omit off-diagonal contributions such as 
(ifJ f2 ' LO ' > )235'32 '°' ) , because they contribute less than 2% of the total. We observe that each 



? {gg^QQ) 



Qo-dependent factor, [in \*f) / In (x)]"""' 3 > ^ s ^ nus suppressed, relative to linear behav- 
ior in a log plot at low Q , for all quasi-octet exponents in the partonic cross section fl58|) . 



This explains why the dashed curve in Fig. 11(a) increases slowly at low Qo compared to 



the dashed curves in Fig. 6(b) 



In Fig. [TU the dashed curve for a pure gluon-induced process in (a) shows a smaller 
gap fraction than the dashed curve for a mixed quark-gluon process in (b) for given Q . 
As we discussed in the previous section, the gap fraction for qq — > QQ rises rapidly at low 
Qo. The term, fg^®®, even with the small assumed probability — q.2 thus explains 

the difference of the dashed curves in Fig. [TU We note that the shape of the solid curves 
for gg — ► Z' — ► QQ in Fig. [TT] corresponds to the largest amount of radiation among the 
resonance processes we have studied in this paper. This is due to the large quasi-singlet 
exponent e^ 9 ^®®^ compared to other exponents. 

The above results, depending on partonic channels and the gauge content of resonances, 
can be understood intuitively in terms of color dipole configurations. In Fig. [TU the color 
dipole structures and the preferred directions of radiation are represented as solid curves and 
arrows, respectively. Each figure describes the surplus of radiation inside a dipole [58-60]. For 



examples, the dipole configurations in Figs. 14(a) and 14(c) allow more radiation into the QQ 
region, which overlaps the gap region, while the configuration for qq — ► G — ► QQ in Fig. 14(d) 
describes enhanced radiation in the qQ and the qQ regions that are outside the gap region. In 



Fig. 14(b) , the color- (anti) symmetric couplings control the dipole configuration, which results 
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(c) qq -» Z' -» QQ 



b 



2, 



(d) qq^G^QQ 



b 



Figure 14: Color dipole configurations for resonance processes, /i(pi) + fiipz) — ► V"' — > 
Q(p a ) + Q{Pb)- Arrows in each figure represent the enhanced direction of radiation. 

in radiation into both the gQ and the QQ regions. We have also predicted the least (most) 
amount of radiation into the gap region for qq — > G — > (<7<7 — >• ^' — »• QQ). These color 
dipole configurations help explain our results for the patterns of radiation, associated with 
the SU(3) color representation of resonances and the partonic channels involved. 

At the end of Sec. El we argued that for a narrow resonance, generally a gauge singlet, 
there is an additional contribution when the radiation energy scale is larger than yTM. In 
this limit, we expect that the color dipole configuration is further dominant since the process 
becomes a pure color singlet decay. This contribution is expected to be absent for an octet 
resonance since the width is in general large. We point out that this may result in more 
radiation for a very narrow singlet resonance process compared to our results in Sees. [5] and 
El leading an even more distinguishable difference from an octet resonance process. The 
details of this case are left for future work. 

7 Conclusion 

In this paper, we have shown that it may be possible to determine the color SU(3) represen- 
tation of resonances from new physics signals by analyzing the distribution of soft radiation 



33 



into a rapidity gap. The results we have found are based on perturbative calculations of 
factorized partonic cross sections, and apply for ordered scales, Aqcd < Q c < \/TM ~ p T 
with energy flow, Q c , resonance mass, M, and decay width, T. These conditions correspond 
to a broad width resonance process. The gap fractions for octet and singlet resonances, 
which we defined in Sec. 12. 3[ show different radiation patterns of energy flow. To obtain 
an analytical form of soft gluon emission in heavy quark pair production, the massive soft 
anomalous dimension matrix for rapidity gap events was introduced. The results, in gen- 
eral, describe more radiation for singlet than for octet resonances, which can be explained in 
terms of color dipole configurations. Especially, for spin-1 resonance production, involving 
only one partonic sub-process, we obtained a relatively large difference in radiation into the 
gap region between color singlet and octet resonances. The results at leading order do not 
require a convolution with PDFs, and are the same at the Tevatron and the LHC in this 
case. At the end of Sec. El we discussed the consequences of a narrow width resonance, 
expecting more radiation for a singlet resonance process compared to our results based on a 
broad resonance assumption. 

In principle, following Ref. [11], given any fixed set of final-state partons, we may calculate 
the distribution of energy flow into any fixed region of rapidity and azimuthal angle rather 
than the simple rapidity gap we defined in this paper. This may allow not only to illuminate 
more distinguishable features of color flow, depending on the gauge content of resonances, 
but also to shed valuable light on the dynamics of QCD radiation itself. 
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Appendix 



In general, we represent the four momentum of a particle of mass mq as 

p^ = (ttit cosh 77, p x , p y , ttit sinh rj) , (A.l) 
where r] is rapidity and where is a transverse mass defined by 

2 2,2,2 2,2 / A 0\ 

m T = mQ + Pa . + Py = mQ + Pr . (A.2) 

In the rj-(f> coordinate system of rapidity and azimuthal angle, four momentum p M can be 
written in terms of quark mass tuq and as 

p M = (y\Jpr J r m 2 Q coshr/, prsin0, pxcos(j), \J Pr + rriQ sinh ij) . (A. 3) 

To describe the process in Eq. (JT]) with the gap geometry, shown in Fig. [1] at lowest order, 
we need five external partons of momenta 



Pi = ^(1,0,0,1) 



p'i = ^(1,0,0,-1), 



A V\ a / „2 , __2 /^ Al A 



Pa = ( \/Pt + m 2 Q cosh [ — 1 , 0, pr, ^/pf + m|sinh I — I 



A ?A n _ /_2 ■ --2 -x. ( A V\ 



Pa = \^/p 2 T + m 2 Q cosh \^—J , 0,-pr, -yp| + m|sinh J 
/c M = /ct (cosh 1/, sin0, cos0, sinh?/) , (A. 4) 

where we are working in the partonic center of mass frame, and where we fix azimuthal 
angles for a quark pair at <p a = and O = ir. Here 7/ is rapidity of soft gluon k. 
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